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Macroscopic quantum tunneling (MQT) is the fundamental dynamics underlying various fields in Physics such 
as the cosmological phase transitions in the very early Universe Q the macroscopic nucleation in ^He — ^ He 
liquid mixtures and phase transitions in nuclear physics and in statistical mechanics. MQT process necessarily 
accompanies many degrees of freedom and is therefore dissipative. This dissipative MQT process is mainly analyzed 
so far in the imaginary time formalism utilizing the instanton method Q . 

However, this imaginary time formalism is not convenient when we study the back reaction of particle creation on 
the tunneling rate. During the system tunnels, it necessarily activates the environmental degrees of freedom which 
couple to the system. The energy of the system is dissipated in this process and the tunneling rate is thought to be 
reduced. If we force to use the imaginary time formalism, we encounter non-Unitary evolution equation for the wave 
function which cannot be normalizable. Furthermore in the imaginary time formalism in general, qualitatively 
different effects such as dissipation and fluctuation cannot be separated in the tunneling dynamics. 

Therefore the real time formalism for the quantum tunneling process is urgently necessary. We explore in this 
paper such possibility utilizing the quantum potential method In this method, the quantum mechanics is exactly 
projected on the classical mechanics with quantum modification of the potential term. This complete one-to-one 
correspondence between the quantum mechanics and the classical mechanics is the starting point of our real time 
analysis. The small probability of quantum tunneling, in this quantum potential method, is not associated with the 
instanton or the saddle point in imaginary time but with a slow-rolling particle on the modified potential in real time. 
Therefore the dynamics of particle creation can be analyzed with the Unitary evolution of wave functions. 

Now we explain the quantum potential method, which was first proposed by Bohm He attempted to construct 
a new interpretation of quantum mechanics(ontological basis) and applied it to the measurement theory of quantum 
mechanics. In this paper we skip the conceptual aspects of this method and directly apply it to the perturbation 
calculation for the modified tunneling rate. 

We start from the one-dimensional quantum system which is described by the Schrodinger equation 
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The polar decomposition of the wave function 5" = i?exp(iT4^/?i) ,with R and W being real functions of a; and i, yields 
an another form of the Schrodinger equation: 

^ + V + Vq = -W (2) 
Am 

and 

2WR + W"R + 2W'R' = 0, (3) 



where Vq is the quantum potential defined by 



^ R" 



The over dot and the prime respectively denote the time and spatial derivatives. The latter equation Eq.(^) simply 
represents the conservation div J(i, x) = of the probability flow J(i, x) defined by 

J(t, x) = {R{t, x)^,Wit, x)'R{t, xf/M). (5) 

The former equation Eq.(||) is the Hamilton- Jacobi equation in classical mechanics if we regard W{t, x) as the action 
function and the total potential Vtot as Vtot = V + Vq- Thus there exists a corresponding classical system for any 
quantum system and this correspondence is exact; it is not any approximation such as WKB. This corresponding 
classical mechanics (effective classical system) is generally different from the underlying classical system (original 
classical system) which is supposed to be quantized to yield the quantum mechanics described by Eq.(^. The whole 
quantum nature is now concentrated on the quantum potential Vq . This correspondence can be extended to systems 
of many degrees of freedom and to relativistic systems. 

This quantum potential method is particularly useful for stationary states. In this case, the variable R becomes 
independent of time, and W — —E. Therefore Vq becomes a time-independent potential. Moreover E — Vtot = 
W'"^ / {2M) > and therefore there appear no classically forbidden region in the configuration space of the effective 
classical system. This fact is particularly interesting when the method is applied to the quantum tunneling process in 
which a particle in the original classical mechanics cannot go over the tunneling region in the configuration space. On 
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the other hand, the quantum fluctuations modify the potential so that a particle in the effective classical system can 
go over the tunneling region. In the effective classical system, the tunneling process is described by a slowly rolling 
particle on the modified potential at the tunneling region in real time. This real time formalism should be compared 
with the ordinary imaginary time formalism using the saddle point method in the path integral (instanton method). 

Now we apply the quantum potential method to the quantum tunneling dynamics. We consider a stationary 
problem for using the advantage of the quantum potential method. First we study a simple one-dimensional model 
with a rectangular potential, which is given by 

{0 for X < (region I) 
Vo for < X < a (region II) , (6) 
for a < X (region III) 

with Vo being a positive constant. The tunneling wave function in each region is given by 

(j)i{x) = Ae''''' + Be-'''^, 0//(a;) = Fe"/^^ + Ge''^, 0///(a;) = Ce*'=^ (7) 



with k = \/2ME/h, (3 — ■\/2M(Vo — E)/h. Smoothness conditions of the wave fmiction at the edges of the potential 
(0 and a) yield 

fF\_ Ce'^" f A_e'3° \ f A\ _ -Ce''"^ / A^e^° - A^e"^" \ 
[gJ^ 2/3 l^A+e-^'^j' [bJ^ Aikfi 1^ ^_A+(e-'3'^ - e^°) j ' 

where A± = /3 ± ik. This wave function yields the quantum potential in each region as is shown in Figs, la and lb. 

Figs, la, lb 

Especially the total potential in region II is given by 



2M [(F-|-/32)cosh(2/3(a-a;)) + (/J^ - fc2)]2 



The exponentially small tunneling probability when /3a ^ 1 is reflected in the fact that in region II, the kinetic energy 
E — Vtot of the classical particle is exponentially small. Actually it is possible to calculate the classical "rolling time" 
over the tunneling region II: 

/■ , r dx 
irou = J at ^ J — 

On the other hand the tunneling probability is given by 

I^P (fc2+/32)2cosh2(/3a)-(/32-/c2)2- ^ > 

When the potential barrier is high [j3a ^ 1), the tunneling probability is proportional to the inverse of the rolling 
time. 

T^-^fE 

P ■ iron « 7= ■ 12 

4VoVVo - E 

Therefore the rolling time can be an another measure of the quantum tunneling. Actually this rolling time is known 
as the "traversal time" . 

For a general form of one-dimensional potential, we can use the WKB approximation with some modification near 
the turning points. We use the connection formula at the turning point a;o where the potential V{x) increases toward 
right {V'{xo) > 0): 
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^ 5 (- \p\dx) 

^ . ,^ 1 
VP 



I E <V side 

(13) 

i= exp ±iij^l:'' pdx+^) mE>V side. 
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where p = ^J2M{E — V{x))/fL. Similar connection form at the turning point with V'{xq) < holds if we replace 
IxQ ^ Ix° ■ important that we cannot drop the exponentially small term which is pure imaginary. This small but 
imaginary term yields non-vanishing kinetic energy of the rolling particle. The graph of the total potential is shown 
in Fig. 2. The total potential in the region near the turning points are estimated by using the airy functions which 
appear when we linearize the potential there. 



Fig. 2 



Before we calculate the precise particle creation effect during the quantum tunneling of the system, we clarify the 
general back reaction effect. The Hamiltonian of the total system we consider is 



H = Hs{x)+H,{x,{y,,}), 



Hs{x) = ^ + Vix), H,ix, {y,,}) = J2 



2to„ 



— '^nOVn + CX 
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We introduce the polar decomposition of the total wave function $ which is assumed to have the factorized form, 

'^{X, {yn}) = 4>s{x)(j)e{x, {y„}), 

0,(x) =i^,(a;)e^^=(-)/^ 
0e(a;, {?/„}) =Y\_Rn{x,yn)expC^iWn{x,yn)/fi)- (15) 

n n 

Putting this decomposition into the total stationary Schrodinger equation for $(a:,{7/n})j we obtain 



{w.s'ix) + j:„w;,) 



2M 
2m„ 



Vix) 



2M RsUnRn 

-2 pa 



+ 2™"'^"02/n + cxyi 



h'' dy^ R„ 
2m„ Rn 



E. 



(16) 



where the prime denotes the derivative with respect to x. 

In order to study the back reaction of the particle creation of the environment on the quantum tunneling of the 
system, we use the quantum potential method as a perturbation from the pure tunneling without particle creation. 
We suppose that the tunneling dynamics of the system alone without back reaction is already known with the system 
wave function i?° exp{iW^ /Ti). Then there is an associated classical motion x{t) derived from Eq.(^). We turn on the 
coupling between the tunneling system and the environment. If this tunneling degrees of freedom is macroscopic and 
semi-classical, then we can define "time" for the environment by using this solution: d/dt = {W^ {x)/M)d/dx. Then 
the wave function of the environment obeys the time-dependent Schrodinger equation with respect to this time 



= He{x{t),{yn})(j)e 



(17) 



in the order 0{h,). A non-trivial time dependence enters from the tunneling solution x{t) which causes particle 
creation. 

Aided by Eq.(|7|), the HJ equation Eq. (|l6|) becomes simpler. 
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(18) 



We suppose that the environment state is well approximated by the Gaussian state. This is guaranteed by the 
Gaussian initial state and the semi- free evolution of the environment. A general Gaussian state is expressed by the 
wave function of the form 



(19) 
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where the real variables and /3„(t) obey the equations of motion: 



, Pn ^ mnUJn{t) . (20) 

m„ nin nin 



Time dependent frequency for the mode n, ujnit) = v^w^q + 2c5;(t)/m„ comes from the motion of the tunneling 
solution X. 

We are interested in the effective motion of the system (x) and therefor should integrate out the environmental 
degrees of freedom. This integration is simply a quantum averaging: 

i^liyn}]) = /n»^^"'^-';f[^>] . (21) 



After this average, the HJ equation finally becomes 



n 

There are two classes of back reaction terms. The first class terms, the second line of Eq. (|2^) , are positive definite 
and they contribute to increase the classical potential. Therefore they systematically reduce the tunneling rate of 
the system. The second class term, the last term in the first line of Eq.(|2^), does not have definite signature and 
may increase or reduce the tunneling rate. It will be interesting that this term is linear in the momentum W'^ which 
reminds us the form of friction. 

We now concretely calculate the particle creation and its back reaction on the tunneling rate of the system. First 
we study the model of rectangular potential. Though the analytic form of the tunneling solution is known from 
Eq.(^, it would be unnecessarily complicated and less general to calculate the normal modes for the time dependent 
frequencies. Therefore we approximate the classical solution a;(t) in the hyperbolic-tangent form: x{t) = a(l-|-tanh(pt)) 
with p = hk/ (oM). This form most faithfully represents the dynamics near the right turning point where the particle 
creation is thought to be maximum^. We would like to solve the evolution equation Eq. (^0|) with this time dependent 
background. Changing the variable from q:„,/3„ to ^„ by ln(^„(t)) = (/3„ + iQ;^J/m„?i, Eq.(|20|) reduces to the form 
^7i{t) + uJnit)'^ ^7i{t) = 0. The solution of this equation can be represented by the hyper-geometric functions. Especially 
the solution which asymptotically approaches to the vacuum for t — )■ — cx3 is given by 

^n{t) = {2uJnQy^^^ exp{iu!n+t + iujn- \n{2 cosh{pt)) / p) X 

X 2Fi[l , ;-(l + tanh(pt))], (23) 

p p p 2 



where LOn± = (^/w^q -I- (4ca)/m„ ± a;„o)/2. Then the typical back reaction terms in Eq.(|2^) are all calculated from 
this solution ^n(0- 



2 / •■ X 2 
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4^.^"-".. . <: ^ . I ^ (24) 

< x5(lne„) a« \2x%\Yiin) 



Since the total back reaction is simply a superposition of each single mode, we concentrate on the back reaction effect 
from a single mode in the following. The result of the numerical calculation for a single mode n is given in Fig. 3. 
Analytically the special case of p = ijJ„o can be easily calculated in the series of e = 2ca/ {mp'^), 



e 



QiU=a = -0.272029- + 0(e2), Q2U=a = 0.14538— + ©(e^). (25) 
a 2a^ 



** This type of bold approximation is inevitable for analytical calculations. 
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Now we evaluate the modified tunneling rate due to the particle creation of a single mode. The total tunneling rate 
can be derived by adding all the single mode effects. The effective Hamiltonian for the system is read from Eq.(|2^) 
excluding the quantum potential, 

^^^^ = 2M + + 4m(^-^°) + W + W 
where po — W^os ■ Then the corresponding Hamilton equation of motion can be brought into a simple form 

Mi + Vlff ^ 0, with 

Veff = V{x] H / poixjax. (27) 



Both classes of back reaction mentioned just after Eq.(|22|) can be expressed in terms of the effective potential Veff 
which is modified from the original potential V{x). In the present case of particle creation from the vacuum, we 
numerically find Qi < as well as Q'l < (Fig. 3). 



Fig. 3 



Therefore we conclude Veff > V{x) and the modification of the potential systematically reduces the tunneling rate. 
This is consistent with the fact that both classes reduce the momentum of the classical particle and increase the 
tunneling time defined in Eq.(^0|). 

The systematic reduction of the tunneling rate comes from the potential shift AV — Veff — V. 



p _ p/'iio particlcN 
^ ~ ^ Vcrcation / ^ 



1 2 . 2MIW 



, 2MaAV, 
exp( ^r^)- (28) 



In the similar way, but using WKB approximation as well as the modification of it near the turning points, 
we can derive the modified tunneling rate for a general potential of the form in Fig 2. If we adopt the same 
strategy to approximate the classical motion by the hyperbolic-tangent form, the parameter p should be given by 
p = (35/^r(2/3)/2r(l/3))(?i/3i/3/Ma), where /3 = -V'{a) with x = a being the right turning point. 

We now summarize our work. We first introduced the quantum potential method in which the quantum mechanics 
is exactly mapped to some classical mechanics (effective classical system) which is generally different from the original 
classical system. When this method is applied to the quantum tunneling process, the effective classical particle can 
slowly roll down the tunneling region in the configuration space. This classical motion in real time corresponds to the 
instanton in the imaginary time formalism. This rolling motion of the tunneling system induces the particle creation 
of the environmental degrees of freedom which were initially in the ground state. As a back reaction of this dissipation 
of energy, the effective potential is increased and therefore the tunneling rate of the system is reduced. 

Several comments on our future work are in order. 

This reduction of the tunneling rate is associated with the initial ground state of the environment. It would be 
interesting to study the particle creation effect for the initial excited state of the environment. In this case, the induced 
particle creation effect would further reduce the tunneling rate while the energy transfer from the environment to the 
system would reduce the rate. 

We started from the classical solution x{t) of the effective classical system, and calculated the modification of the 
effective potential. Finding a new solution for this effective potential, we can repeat this quantum potential method. 
Then we may be able to find a series of solution which asymptotically approaches the true wave function including 
the full back reaction of particle creation. This iteration method would be another approximation paradigm which is 
free from any turning points since there is no classically forbidden region in the effective classical system. 

We can directly apply our method to the quantum cosmology. Especially our method would reveal back reaction of 
the cosmological particle creation on the tunneling solution of the Wheeler-DeWitt equation. We will soon reported 
these work. 
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Figure captions 

Fig. la: The rectangular potential V (dashed graph), and the associated total potential Vtot in the tunneling region 
(sohd graph). We set the energy E = 2 (dotted graph), a = 1, Vq = 4 in the units h = 1,M = 1. 

Fig. lb: Same as Fig. la but in global view. 

Fig. 2: The quadratic potential V{x) = 1 — 8x{x — 1) (dashed graph) and the associated total potential Vtot (solid 
graph) calculated from WKB approximation with modification by Airy function near the turning points. We 
set the energy E = 1 (dotted graph) in the units h = 1,M = 1. 

Fig. 3: The particle creation effect on the tunneling rate is fully represented by the effective potential Vctf{x) (solid 
graph) which is calculated for the rectangular potential model. We set the parameter as£' = 2,a=l,Vb = 
4:,m = l,u)Q = l,c = 0.15 in the units U = 1, M = 1. 
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